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Abstract 

We formulate a generalized mean-field theory of a mixture of fermionic and bosonic atoms, 
in which the fermion-boson interaction can be controlled by a Feshbach resonance. The theory 
correctly accounts for molecular binding energies of the molecules in the two-body limit, in con- 
trast to the most straightforward mean-field theory. Using this theory, we discuss the equilibrium 
properties of fermionic molecules created from atom pairs in the gas. We also address the forma- 
tion of molecules when the magnetic field is ramped across the resonance, and present a simple 
Landau-Zener result for this process. 

PACS numbers: 
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I. INTRODUCTION 



The use of magnetic Feshbach resonances to manipulate the interactions in ultracold 
quantum gases has greatly enriched the study of many-body physics. Notable examples 
include the crossover between BCS (Bardeen-Cooper-Schrieffer and BEC (Bose-Einstein 
Condensate [2j, yf ) superfluidity in ultracold Fermi gases [4 a ]6J , and the "Bose Nova" 

n rf n nnnn 

collapse in Bose gases [7|. Recent experimental developments, |8j, |9|, [lGj, [III, ll2l, [13j have 
enabled the creation of an ultracold mixture of bosons and fermions, where an interspecies 
Feshbach resonance may introduce a rich source of new phenomena. From the theoretical 
point of view, studies of Bose-Fermi mixtures to date have been mostly limited to non 



resonant physics, focu sing mainly on mean field effects in trapped systems 
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24j . phases in optical lattices 
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13, 



|30j,or equilibrium 



studies of homogeneous gases, focusing mainly on phonon induced superfluidity or beyond- 

S, [37! ■ Pioneering theoretical work on the resonant gas 



mean-field effects 
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include Ref. 38j, in which a mean-field equilibrium study of the gas is supplemented with 
a beyond-mean-field analysis of the bosonic depletion; and Ref. [39], where an equilibrium 
theory is developed using a separable-potentials model. 

The aim of this article is to develop and solve a mean-field theory describing an ultracold 
atomic Bose-Fermi mixture in the presence of an interspecies Feshbach resonance. This 
goal appears innocuous enough at first sight, since mean-field theories for resonant Bose- 
Bose 40, |4lJ and Fermi- Fermi 42|, |43|, |44j gases exist, and have been studied extensively. 
In both of these theories, the mean-field approximation consists of considering the bosonic 
Feshbach molecules as being fully condensed, and this greatly simplifies the treatment, since 



45|. 



the Hamiltonian reduces to a standard Bogoliubov-like integrable form 

The fundamental difference between these examples and the Bose-Fermi mixtures is 
that in the latter the Feshbach molecules are fermions, and therefore their center of mass- 
momentum must be included explicitly. The most obvious mean-field approach consists in 
considering the atomic Bose gas to be fully condensed. However, as we will show below, res- 
onant molecules are really composed of two bound atoms, which spend their time together 
vibrating around their center of mass. It follows that outright omission of the bosonic fluc- 
tuations of the atoms, disallows the bosonic constituents to oscillate (i.e. fluctuate) at all, 
and therefore this leads to an improper description of the physics of atom pairs. 
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This article is organized as follows. In section [TT] we introduce the field theory model used 
to study Feshbach resonances, describing briefly the parametrization used, and outlining the 
exact solution of this model in the two-body limit. The section ends with a test of this two- 
body theory, by comparing the binding and resonance energies predicted by the model and 
the virtually-exact analogues obtained from two-body close coupling calculations. Section HTT1 
introduces the simplest mean-field many-body theory of the gas, obtained by disregarding all 
bosonic fluctuations. The solution of the theory is outlined, and its limitations highlighted. 
In spite of these limitations, mean-field theory provides a useful language for dealing with 
the problem, the utility of which will persist even beyond the limits of applicability of the 
theory itself. 

Finally in section [IV] we introduce our generalized mean-field theory, which is, in short, 
similar to the mean-field theory described in section HTT| but with the notable improvement of 
using properly renormalized molecules as building blocks, instead of their bare counterparts. 
This approach is not trivially described in the Hamiltonian formalism, where substituting 
dressed molecules for free ones would lead to double counting of diagrams. In this section 
we therefore shift to the Green-function/path-integral language, where this double-counting 
can be avoided quite easily. Finally we proceed to the numerical solution of this theory, and 
note that for narrow resonances the results are consistent with their mean-field equivalents. 
This encourages us to develop a simple theory to study the molecular formation via magnetic 
field ramps, and, using an approach based on the Landau-Zener formalism Q, I^tJ , we derive 
analytic expressions in section IVl 

Throughout this article, we work with zero temperature gases, in the free space ther- 
modynamic limit. These are limitations which render the results obtained here hard to 
directly compare with experimental results. One of the main possible future directions of 
this work should include solving the same problem in a trap, and generalization to higher 
temperatures. 

II. THE MODEL 

We are interested primarily in the effects of resonant behavior on the otherwise reasonably 
understood properties of the system. To this end we use a model which has become standard 
in the last few years. This model has been useful in studying the effect of resonant scattering 
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in Bose 
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49(1 and Fermi 
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mixture, this model has been used in Ref. 
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,and 



55J gases. In the case of the bose-fermi 



56 



571 ] . We refer to these works for 



further details about the origin and justification of the Hamiltonian we use here, and for 
details on the solution in the two body regime. 
The Hamiltonian for the system reads: 



H = H + Hj, 



(1) 
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Here the annihilator operators for, respectively, fermions and bosons, c p is the 



annihilator operator for the molecular field 
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551 ]; 7 = Airab/mi, is the interaction 
term for bosons, where a& is the boson-boson scattering length; Vb g , u, and g are parameters 
related to the Bose Fermi interaction, yet to be determined; the single particle energies are 
e a = p 2 /2m a , where m a indicates the mass of bosons, fermions, or pairs; and V is the volume 
of a quantization box with periodic boundary conditions. 

The first step is to find the values for V& s , u, g, in terms of measurable parameters. We 
will, for this purpose, calculate the 2-body T-matrix resulting from the Hamiltonian in Eq. 
[2] Integrating the molecular field out of the real time path integral, leads to the following 
Bose-Fermi interaction Hamiltonian: 



H 



2body 



V [ Vb9+ E 



v 



(3) 



This expression is represented in center of mass coordinates, and E is the collision energy of 
the system. From the above equation we read trivially the zero energy scattering amplitude 
in the saddle point approximation: 



(Vi 



bg 



), 



(4) 



which corresponds to the Born approximation (This is akin to identifying the scattering 
amplitude / = a sc in the Gross-Pitaevskii equation, where the interaction term would be 
^j a b )• We emphasize that this approximation is only valid in the zero energy limit, and it 
does not, therefore, describe the correct binding energy as a function of detuning. However, 
with this approach we obtain an adequate description of the behavior of scattering length as 
a function of detuning, which allows us to relate the parameters of our theory to experimental 
observables via the conventional parametrization [|5(3] 

T ^ = ( 1 " (5^o) • < 5 > 

where a\, g is the value of the scattering length far from resonance, A# is the width, in 
magnetic field, of the resonance, m&/ is the reduced mass, and Bq is the field at which the 
resonance is centered. 

The identification of parameters between Eqns., (HI) and (j5J) proceeds as follows: far from 
resonance, \B — B \ >> A#, the interaction is defined by a background scattering length, 
via VL = 27ra6g . To relate magnetic field dependent quantity B — B to its energy dependent 
analog u, requires defining a parameter 5b = dv/dB, which may be thought of as a kind 



of magnetic moment for the molecules. It is worth noting that v does not represent the 
position of the resonance nor the binding energy of the molecules, and that, in general 5 b 
is a field-dependent quantity, since the thresholds move quadratically with field, because 
of nonlinear corrections to the Zeeman effect. For current purposes we identify 5b by its 
behavior far from resonance, where it is approximately constant. Careful calculations of 
scattering properties using the model in Eq. (j2J), however, leads to the correct Breit-Wigner 
behavior of the 2-body T- matrix j5s|. 

Finally we get the following identifications: 

T/ 27ra bg 

g = ^jV bg 5 B &B 
u = 5 B (B-B ). 

(6) 

I — I 

Synthesizing the approach descibed in [56(], and diagrammatically represented in Fig. (H 
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we can obtain the exact two-body T-matrix of the system by solving the Dyson equation 



T = gD°g + gD°g U gD°g + gD°g U gD°g U gD°g + ■ 
= gDg 

where T is the T-matrix for the collision, and which has formal solution 

,2 



(7) 



T = gDg 
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(D*)- i - g m 

These quantities take the explicit form 

du dp 1 



(8) 



11(E) 



2vr (2vr)3 + W 

1 



X 



(9) 



(E-tuj-^+iOr 



m 



3/2 
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m bf A 



(10) 



V2n V 

where my is the boson-fermion reduced mass, and A is an ultraviolet momentum cutoff 
needed to hide the unphysical nature of the contact interactions. Note that D° represents 
an effective molecular field, accounting for the fermion-boson background interaction, and, 



as described in detail in 



561 ]. it is obtained by integrating the original molecular field, 



and performing a Hubbard-Stratonovich transformation 59| to eliminate the direct boson- 
fermion interaction in favour of the effective molecular field D°. 



Toe 

f = 



+ 



+ 



g D° g + g D°gU g D°g 



g D° g II g D° g ilg D° g 



+ 



+ 



g D g 



FIG. 1: Feynman diagrams representing the resonant collision of a fermion and a boson. Solid 
lines represent fermions, dashed lines bosons, and double solid-dashed lines represent the effective 
composite fermions. 



6 



Regularization of the theory [42j, |56| is obtained by the substitutions 




(11) 



Finally the two-body T-matrix takes the form 




-i 



T(E) 



1 



(12) 



Vbg + 



A. Poles of the T-Matrix: Testing the Model 

Bound states and resonances of the two-body system are identified in the structure of 
poles of the T-matrix (Eq fT2|) , This is illustrated in Figure |2j where real and imaginary parts 
of the poles' energies are plotted as a function of magnetic field. The resonance portrayed 
in the figure is the 544. 7G resonance present in the 1 9/2, — 9/2) |1, 1) state of 40 K- 87 Rb. For 
B < 544. 7G (corresponding to detunings v < 0) the two-body system possesses a true bound 
state, whose binding energy is denoted by the solid line. In this case, the pole occurs for 
real energies. This bound state vanishes as the detuning goes to zero, where the resonance 
occurs. 

For positive detunings, v > 0, on the other hand, the poles are complex, and the inverse 
of the imaginary part is proportional to the lifetime of the metastable resonant state. In 
this regime, there is no longer a true bound state, but there may be a scattering resonance, 
indicated in Fig. [2] by a thick dashed line. This resonance appears for magnetic fields 
B > 544.7 for this particular resonance, well before the disappearance of the bound state. 
This value is highly dependent upon the value of the background potential. We will see in 
section HVl that for V& s = 0, the resonance actually appears at positive detunings. In the 
case of 40 K- 87 Rb, Vb g < 0, implying that there is a weak potential resonance in the open 
channel which interferes with the closed- channel resonance, and causes it to cross the axis 
at negative detunings. For Vbg > ([60]) the positive background scattering length is set 
by a bound state in the open channel, which does not affect the resonance states, but which 
interferes with the bound state at negative detunings. 
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FIG. 2: The top panel shows the scattering length versus magnetic field for the 544.7 G resonance, 
present in the 1 9/2 , — 9/2) |1, 1) states of the 40 K- 87 Rb collision. The bottom panel shows the poles 
of the model two-body T-matrix (Eq. fl~2j) parameterized for the same resonance, as a function of 
magnetic field. Thick solid and dashed lines denote the real parts of relevant poles, representing 
bound and resonance states respectively. The thin dashed lines are the real parts of unphysical 
poles. The empty circles and squares, represent the position of the resonance and the bound state, 
obtained via a virtually exact close coupling calculation, and are presented to show the level of 
accuracy of the model. 
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The thin dashed lines in Fig. [2] are physically meaningless solutions to the Schrodinger 
equation, in which the amplitude in the resonant state would grow exponentially in time, 
rather than decay. These poles do not therefore identify any particular features in the energy- 
dependent cross section of the atoms, and will not modify the physics of the system. Finally 
Fig. [2] contains data obtained from virtually exact close-coupling calculations which show 
the extent of validity of the model. For the purposes at hand this agreement is sufficient. 

It should be noted that the agreement is not as good for positive background scattering 
length systems, since the open-channel bound state determining this scattering length is not 
adequately described by the model, which treats the background physics as an essentially 
zero range interaction. This implies that the relation between the background scattering 
length and open channel bound-state energy is exactly Ej, = l/2fial g , while in the physical 
system this relation depends on the details of the interaction potential. This problem has 
been addressed in the literature [6l| . but no treatable field theory has yet been proposed. 



B (G) A B (G) 8 B {K/G) a bg {a.u.) 



492.49 0.134 3.624 x 10~ 5 -176.5 

544.7 3.13 1.576 x 10~ 4 -176.5 

659.2 1.0 2.017 x 10~ 4 -176.5 



TABLE I: Parametrization of the three main Feshbach resonances used in this thesis. All Three 
resonances are in the 1 9/2, — 9/2) |1, 1) states of the 40 K — 87 Rb collision 



III. MEAN-FIELD THEORY 

In this section we introduce the many-body physics of the system, by first analyzing it in 
a mean-field approach. Because of the statistical properties of the system, we will see right 
away that mean field theory does not recover the correct two body physics in the low density 
limits. In spite of this substantial weakness, however, the approach has several qualitative 
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features which persist even in the improved theory that we introduce below. Furthermore, 
since the model is exactly solvable, it will allow us to develop a language which will help 
us to understand the problem in simpler terms, and to identify some small physical effects, 
which, when ignored, can greatly simplify the beyond mean-field approach presented in the 
next section. 

A. The Formalism 

Starting with the Hamiltonian described by equation ([I]), we obtain the mean- field Hamil- 
tonian by substituting the boson annihilator b by its expectation value = (b), a complex 
number. The number operator tf p b p therefore becomes \<fr\ 2 = iV&, where N b is the number of 
condensed bosons. The grand canonical Hamiltonian therefore becomes 

H = E b +} j (ep - fj,f + VbgUb) a\a p (13) 
p 

V 

v 

where rib is the density of condensed bosons, Eb/V = •ynl—fibnb is the energy per unit volume 
of the (free) condensed bosons, a constant contribution to the total energy of the system, and 
A*(6,/,m) are the chemical potentials. These are Lagrange multipliers that serve to keep the 
densities constant as we minimize the energy to find the ground state. In the following we will 
drop the volume term, absorbing it in the definition of the creator /annihilator operators, 
such that the expected value of the number operator represents a density, instead of a 
number. 

Before proceeding with the analysis of this Hamiltonian, we introduce the set of self- 
consistent equations we wish to solve. To this end we define the quantities representing 
the total density of bosons (fermions) in the system, at detuning v — ► oo. At finite detunings 
some of these atoms will combine into molecules, and the densities will be denoted as nibj, m ) 
for bosons, fermions and molecules respectively. 

The system, therefore, is described by six quantities, namely three densities and three 
chemical potentials, which require six equations to determine. These equations, which can 
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be derived by number-conservation constraints and energy minimization arguments, are: 



Tif + n m — n°f — 


(14.a) 


n b + n m - n° b = 


(14.b) 
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(14.c) 


d n 
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dfji m 


(14.d) 


dVL 
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(14.e) 




(14.f) 



where Q = (H)/V is the Gibbs free energy. Equations ( 11 4. aft and (114. bl) follow from the 



simple counting argument that for every molecule created, there is one less free boson and 
one less free fermion in the gas. Equations I14.cl and ll4.dl are simply the Lagrange multiplier 
constraint equations, equation (114. el) follows from the mean-field approximation, whereby 
the bosonic field is simply a complex number, and minimization of the energy can therefore 
be done directly. Finally equation ( 114. fl) is the law of mass action, which follows from the 
fact that to make a molecule it takes one free atom of each kind. 

The next step is to write down Q for the system, by taking the expectation value of the 
Hamiltonian in equationf [T3j). obtaining 

= E b /V + J2{e*-Hf + V bg n b ) rj f (p) (15) 
p 

p 

p 

where rjf(p) = (d) p d p ) and T] m (p) = (c^c p ) are the fermionic and molecular momentum dis- 
tributions, T] m f(p) = (Cptip) is an off-diagonal correlation term arising from the interactions 
in the system, and the densities are given by nf jm>m f = J A?7/,m,m/(p)- Equations I14.atfl~4~fl 
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then read 



rif + n m — = (16. a) 

rib + n m - n Q b = (16. b) 

n f,m,mf = I ^Vf,m,mf(p) (16. c) 



g n mf - fi by /n~ b + >yn b /2 = (16. d) 

Hb + Hf = Hm (16-e) 

The remaining task is now to find expressions to calculate the expected values r if^ .mfiv)- 
To this end we follow a Bogoliubov-like approach, similar to that described in 38| . The 
mean-field Hamiltonian is bilinear in all creation/annihilation operators, which means that 
it can be diagonalized via a change of basis, whereby introducing the operators 

a p = A a a p + C a c p (17) 
f3 p = Apcip + Cp&p, 

for some appropriately chosen coefficients A a ,p and -B Ql/ 3, the Hamiltonian will read 

YsMpY^A (18) 

p v 

At this point we note that the Hamiltonian is just a separable sum of free-particle Hamil- 
tonians, where the free particles are fermions, with dispersion relations X a ^(p). We can 
readily write down the distribution 

V*AP) = ©(-Va(p)) (19) 

where 6 is the step function, and calculate the densities n a ^. The step function could 
be replaced by the free Fermi distribution for non-zero temperatures, but the mean-field 
assumption that all bosons are condensed would no longer hold. If these were ordinary 
free fermions with dispersion p 2 /2m — fi, equation f|T9|) would reduce to the standard zero- 
temperature Fermi distribution. We will see below that \ a ,p{p) are dispersion relations of 
quasi-particles that are a mixture of atoms and molecules. 

Below we show how these ideas, together with equation (I14.a[) - Il4.fl give us the tools 



we require to calculate the observable atomic and molecular densities as a function of the 
chemical potentials. 
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To illustrate more explicitly the diagonalization procedure we define the vectors 

A=\ ap ) ^ = ( aj fit ), (20) 



and 

a 



B= yHj st =(*{#). ( 2i ) 

whereby the Hamiltonian can be written as A^HA, and B^H'B, where 



and 



W = | A " W ° I (23) 



o v(p) 

Diagonalizing if, we get the two eigenvalues 



w = Mg) + Mp) (24) 



± 

2 

where we have defined hf(p) = (e£ — /// + Vj, g n&J, and h m {p) = (e^f + v — fM m ), and the 
unitary eigenvector matrix 

U=( AaBa ). (25) 

The transformation in eq. [IT] can then be written as A = U^B, and its inverse B = UA. 

Our goal now is to write the densities r)m,f,mf(p) m terms of the known densities r) a> p(p). 
In component notation, (where Ai = a p , etc.), we can write 

(AlA) = (BjU Jt (U%B k ) = UyUMBfa), (26) 

where we have used the fact that since the Hamiltonian is diagonal in the B basis, then 
(BjB k ) = {B]B 3 )5 3k . 

Using this formalism we obtain the relations: 

77/(p) = \A a \ 2 V aip) + \B a \ 2 r ]l3 i P ) (27) 
Vf(p) = \Ap\ 2 r] a (p) + iBpfrj^p) 
(p) = A* a A pVa (p) + B* a Bp Vp (p). 
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Using these expressions in conjunction with eqs. U6.aH16.el will then allow us to compute the 
equilibrium properties of the system. 

This simplified, mean-field version of the solution can only approximately reproduce the 
energies of atomic and molecular states, as is shown in figure [3J In this example, we have 
assumed a uniform mixture of 40 K and 87 Rb atoms with densities 8.2 x 10 14 cm -3 and 
4.9 x 10 15 cm" 1 , respectively. These densities correspond to the central density of each 
species assuming it is confined to a 100 Hz spherical trap. Far from the resonance, these 
energies asymptote to zero (representing the atomic state) and to the detuning (representing 
the bare molecular state that went into the theory). Near zero detuning, these levels cross 
owing to the coupling term g^/nb in Eqn. ( 1221) . The size of this crossing is therefore larger the 
larger the bosonic density is. In the crossing region the eigenstates do not clearly represent 
either atoms or molecules, but linear combinations of the two. 




-10 -8 -6 -4 -2 2 4 6 8 10 



FIG. 3: The thick solid lines represent the "renormalized" mean-field energy levels X a ^(kf) , while 
the thin dashed lines represent their bare counterparts. 
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B. Mean field: noninteracting case 

To better understand the structure of the mean-field theory, in this section we detail its 
results for a noninteracting gas, by setting g = Vb g = 7 = 0. We contrast two different 
physical regimes, based on the ratio of bosons to that of fermions, = n^/uf. 

In the case of high fermion density, we set = 0.6, and plot chemical potentials and 
populations of the various states in Figure HI Consider what happens in an infinitely slow 
ramp from positive detuning (no molecules) to negative detuning (introducing bound molec- 
ular states). For large positive detuing, the fermionic and molecular chemical potentials are 
the same, since the chemical potential of the condensed bosons vanishes. For large enough 
detuning, the molecular chemical potential remains below the detuning, so it is energeti- 
cally unfavorable to make molecules. When the detuning dips below the chemical potential 
(detuning a in the figure), fermions begin to pair with bosons to make molecules (see popu- 
lations in lower panel). This process continues until all the bosons are consumed (detuning 
b), at which point the populations stabilize. For detunings less that this, there remain both 
fermions and fermionic molecules in the gas, and there are two Fermi surfaces present. Be- 
cause the internal energy of the molecules continues to diminish at lower detuning, so does 
the molecular chemical potential. The two Fermi surfaces therefore split from one another, 
although the relative population of the two fermions is fixed. 

By contrast, the case where the bosons outnumber the fermions is shown in Figure 
where we have set the density ratio to r&/ = 6. As in the previous case, no molecules are 
generated until the detuning drops lower than the chemical potential of the atomic Fermi 
gas. Since there are enough bosons to turn all the fermions into molecules, there there are no 
fermionic atoms at sufficiently negative detuning, and the gas possesses only a single Fermi 
surface. The chemical potential of the remaining bosons is still zero, since these bosons are 
condensed. Formally, then, the chemical potential for atomic fermions is negative, meaning 
that their formation at negative detuning is energetically forbidden. 



15 




FIG. 4: Equilibrium chemical potentials (top panel) and populations (bottom panel) as a function 
of detuning for a non interacting gas with r&/ = .6. The solid lines represent fermions, dashed lines 
molecules, and dashed-dotted lines bosons. The dotted lines in the top panel represent the bare 
molecular and fermionic internal energies, respectively v and 0. The vertical lines labeled a) and 
b) are discussed in the text, and represent the detuning at which molecular formation begins and 
ends, respectively 
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FIG. 5: Equilibrium chemical potentials (top panel) and populations (bottom panel) as a function 
of detuning for a non interacting gas with r^f = 6. The solid lines represent fermions, dashed lines 
molecules, and dashed-dotted lines bosons. The dotted lines in the top panel represent the bare 
molecular and fermionic internal energies, respectively u and 0. 

C. Mean field: interacting case 

In most experimental circumstances, the density of bosons is larger than that of fermions, 
since condensed bosons cluster to the center of the trap, whereas fermions are kept away by 
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Pauli blocking. We therefore focus on this case hereafter, setting the bose and fermi densities 
to rib = 4.9 x 10 15 cm -3 and n/ = 8.2 x 10 14 cm -3 . The coupling term g-Jni in equation fl22|) 
is the perturbative expansion parameter for the problem, and since it has units of energy, 
it must be compared with the characteristic non-perturbed energy of the gas, which in this 
case is Ef. Also since in the perturbative expansions it always appears squared (see eq. Ej), 
we can define the unitless small parameter for the system as e$M — 9 2 nb/Ej = g 2 rif/ Ejr^f, 
For the 492G resonance in table [U we have €sm = 6.35 x 10~ 2 r b f, and since r b f = 6, the 
small parameter is of order 0.1, appropriate for perturbative treatment. 

Figure [6] shows the equilibrium chemical potentials for the system obtained via a self- 
consistent solution of equations ( 1281) and I16.alll6.el This figure is qualitatively similar to 
the corresponding non-interacting result in Fig. [5j but contains important differences. First, 
the nonzero boson-boson interaction generates a nonzero bosonic chemical potential fib that 
breaks the degeneracy between the molecular and fermionic chemical potentials. In physical 
terms, this means that there is an energy cost in maintaining bosons unpaired, and therefore 
we need to take this into account in the kinematic analysis. Namely, to make molecules 
energetically favorable no longer requires a detuning v such that v = iif, but now requires 
v = fif + fif,- The net result is to shift the chemical potential up and to the left by an amount 
fj,b, and to shift the molecular population curve to the left by this amount. 

A second difference is that molecule creation takes place more gradually as a function of 
detuning in the interacting case. This is simply the result of the avoided crossing smearing 
out the molecular energy. 

Finally, we exploit the simplicity of the mean field approach to test some approximations 
that will simplify the beyond-mean-field approach in the next section. These approxima- 
tions have been tested numerically, and they give corrections of the order of .1% or less 
in calculated molecular populations for all regimes of interest here. The approximations 
are: i) incorporate the boson-boson interaction 7A" 3 by shifting the detuning and chemical 
potential as discussed above; ii) disregard the background scattering between bosons and 
fermions, i.e., set Vb g = 0, since this interaction is dominated by its resonance part; and iii) 
disregard the correlation function (rjbf{p)) ((analogous to the boson polarization operator 
in the Green function formalism), since its contribution to fib is much smaller than that 
of 7n 3 . It is difficult to directly verify the validity of these approximations in the beyond- 
mean-field approach. Nevertheless, we expect that these approximations remain valid, since 
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FIG. 6: Mean-field equilibrium chemical potentials (top panel) and populations (bottom panel) 
as a function of detuning. The solid lines represent fermions, dashed lines molecules, and dashed- 
dotted lines bosons. 



the generalized mean field theory is, after all, a mean-field theory at heart. 
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IV. GENERALIZED MEAN-FIELD THEORY 



In section II III we reached the conclusion that the mean-field approach to the resonant 
Bose-Fermi system does not properly account for the correct two body physics of the system. 
In this section we wish to improve on this, by introducing a generalization to mean-field 
theory, via an appropriate renormalization of the molecular propagator, which is able to 
reproduce the correct two-body physics in the low-density limit. To accomplish this, we will 
have to abandon the Hamiltonian treatment of the previous section, in favor of a perturbative 
approach based on the Green's function formalism, much as was done for two bodies in HT1 
Throughout this section we use the approximations made above, namely, ^n b = V bg = 
(Vbf(p)) = o. 

We begin by recasting the rn field result from Sec. IHIl in the language of Green functions. 
The self-consistent Dyson equations that describe this system are 

G™ F (E,P) = G F (E,P)+g 2 n b D°(E,P) G^ IF (E,P) 

D MF (E, P) = D°(E, P) + g 2 n b G° F (E, P) D MF {E, P), (28) 

where the free propagators are simply 

D°- 1 



u _ f AT(p) + ^ S i gn (f (p)) 

qo = \ (29) 

F ^~i F {p) +iV sign(f(p))' 

and i M ' F {p) = (e^' F — Mm,/) These propagators are described diagrammatically in Fig. 
They represent the fact that a free fermion may encounter a condensed boson and associate 
with it, temporarily creating a molecule; or that a free molecule may temporarily split into 
a fermion and a condensed boson. Self-consistency ensures that these processes may be 
repeated coherently an infinite number of times. We neglect the bosonic renormalization 
equation <p MF = + g 2 n b G° F (E, P) D°(E, P), whereby a condensed boson may pick-up a 
fermion to create a molecule; this is equivalent to the condition (rjbfip)) = 0. 
Solutions to these equations take the form 

1 



G»*{E,P) 



G° F (E,P)-i-g 2 n b D°{E,P) 



° MF ^ P) - D^P^-g^G^PY ^ 



20 



— • + 



™ r r 




FIG. 7: Feynman diagrams included in the mean- field theory. Thin (thick) solid lines represent 
free (renormalized) fermions, thin double dashed-solid lines represent free molecules, and thick dou- 
ble dashed lines represent renormalized molecules. The little lightning bolts represent condensed 
bosons, whereby the arrow indicates whether they are taken from or released into the condensate. 

Using the definitions of G F (E,P) and D°(E,P) from Eq. [29J, we can find the poles cor- 
responding to many-body bound states. In this case, it can be shown that the poles are 
exactly the mean-field eigenvalues \ a ,p(p) from section [Til Moreover, the equations (|30l) are 
symmetric with respect to interchange of G F and D, which implies that both renormalized 
green functions have the same poles, and the same residues. We can therefore study the 
properties of the fermions by only looking at the molecules. This is not completely surpris- 
ing, since, given that the condensed bosons are relatively inert, every molecule corresponds 
exactly to a missing fermion, and vice-versa. 

The most important deficiency of this mean-field approach is that it only allows molecules 
to decay into a free fermion and a condensed boson, disregarding the possibility that the 
bosonic byproduct may be noncondensed. We must allow noncondensed bosons somehow, 
and yet these bosons make a perturbation to the result, as seen by the following argument. 
The fundamental mean field assumption is that the gas is at zero temperature, and therefore 
the noncondensed population should be negligible at equilibrium. Furthermore, if a molecule 
is composed of a zero-momentum boson and a fermion from the Fermi sea, dissociating into 
a noncondensed boson implies that the outgoing fermion would have momentum lower than 
the Fermi momentum, an event which Pauli blocking makes quite unlikely. Therefore, if a 
molecule does indeed decay yielding a non-condensed boson, it should immediately recapture 
the boson in a virtual process such as that described in Fig. HJ It is only convenient that these 
events are exactly the kind of events which will correctly renormalize the binding energy of 
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the molecules, leading to a theory which will reproduce the exact two-body resonant physics. 
The Dyson equation describing this generalized mean-field theory are: 

G GMF (E,P) = G° F (E,P) + g 2 n b D(E,P) G GMF (E,P) 
D GMF (E,P) = D(E,P) + g 2 n b G F (E,P) D GMF (E,P). 

(31) 

Here we have replaced the free propagator D° by the renormalized molecular propagator 
D from equation (jSJ). A diagrammatic representation of this theory appears in Fig. [SJ By 
analogy with the mean field version, the solution to these equations are: 

1 



G GMF (E,P) 



G F (E,P)-i-gin b D(E,P) 



DGMF ^ = Dppikw) - (32) 

These equations preserve the symmetrical nature of the mean-field theory described above, 
and also the avoided crossing of atomic and molecular levels. This is demonstrated in Figure 
|9l where we reproduce the P = kf pole of D from Fig. [2] as dashed lines. We also present 
in this figure the corresponding poles for the generalized mean-field theory, as solid lines. 
As in figure [3] we note the splitting in two energy levels, avoiding each-other around v — 0. 
The fundamental difference in this case is that the molecular curve does not asymptote to 
the bare detuning, but rather to the correct molecular binding and resonance energies. 

Studying the equilibrium properties of the system is now a matter solving the self- 
consistent set of equations fll6.aH16.eI) . while setting rj m f and A equal to zero. To do this we 
first need to extract the distributions r)f }Tn from the Green functions D GMF and G GMF . To 
avoid taking a distracting detour here, we refer the reader to appendix |A] for details. As 
in the previous section, we will consider a mixture composed of a free gas of fermionic 40 K 
atoms, with a density of 8.2 x 10 14 cm -3 , and a gas of condensed 87 Rb bosons with density 
4.9 x 10 15 cm -3 (corresponding to the respective Thomas-Fermi densities of 10 6 atoms of 
either species in the center of a 100Hz spherical trap). 

Figure [TD] shows the equilibrium molecular population as a function of detuning, for the 
492. 5G resonance. For the densities assumed, the mean-field parameter €sm — g 2 n b /E"j w 
0.4 is indeed perturbative. For this narrow resonance, the agreement between mean-field and 
generalized mean-field is quite good, and we could as easily have used the bare molecular po- 
sitions to calculate this quantity. However, the situation is completely different for the wide 
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FIG. 8: Feynman diagrams included in the generalized mean- field theory. Like in the mean field 
case (Fig. [7J, thin (thick) solid lines represent free (renormalized) fermions, thin double dashed- 
solid lines represent free molecules, and thick double dashed lines represent renormalized molecules. 
The little lightning bolts represent condensed bosons, whereby the arrow indicates whether they 
are taken from or released into the condensate. The novelty here in the inclusion of the 2-body 
dressed molecules from Fig. [TJ 

resonance at 544. 7G, for which the equilibrium molecular populations are shown in Fig. [TT1 
Here the "perturbative" parameter has the value 6sm = 38.7, and is not perturbative at all. 
For a given small detuning, the simple mean-field approximation would greatly overestimate 
the number of molecules in the gas at equilibrium. The more realistic generalized mean field 
theory accounts for the fact that the actual molecular bound-state energy is higher than 
the bare detuning. This fact in turn hinders molecular formation, according to chemical 
potential arguments analogous to those in Sec. IIII Bl 

V. MOLECULE FORMATION 

Moving beyond equilibrium properties, we are also interested in the prospects for molecule 
creation upon ramping a magnetic field across the resonance. In reference [56], the mean-field 
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FIG. 9: The thick lines represent the poles of D (kf), and correspond to molecular and fermionic 
energies. The dashed lines represent the molecular two body poles, corresponding to the molecular 
binding and resonant energies obtained disregarding the background interaction [571 ]. As in the 
mean-field theory case, the effect of the interaction with the condensate is to create an avoided 
crossing between the atomic and molecular states. The "bulge" in the upper solid curve is a 
consequence of Pauli blocking which has the consequence of favouring molecular stability. See [5^] 
for more details. 



24 




FIG. 10: Equilibrium molecular population as a function of detuning for the narrow 492. 49G 
resonance. The solid line represents results obtained via the generalized mean field theory presented 
in the text, while the dashed-dotted line represents the mean field results. 



equations of motion for the system at hand were derived, as follows: 

d 

d 

h—r] F (p) = -2g Qm((f)r] MF (p)) 
d 

h—r)M[p) = 2g ^sm((j)r] MF (p)) 
d 



ifr-^;VMF(p) = [4 - tp - v + V bg \<p\ 2 ] t}mf{p) ~ 
g<f)* (r] F (p) -t]m(p)), 



(33.a) 
(33.b) 
(33.c) 

(33.d) 



25 




100 



FIG. 11: Equilibrium molecular population as a function of detuning for the wider 544. 7G reso- 
nance. The solid line represents results obtained via the generalized mean field theory presented 
in the text, while the dashed-dotted line represents the mean field results. 

where r] F (p) = (a) p a p ) is the the fermionic distribution, rj M (p) its molecular counterpart, 
and Pm,f = / ^2P 2 Vm,f(p) the fermionic and molecular densities. Similarly !Jmf(p) = (c p a p ) 
and is the distribution for molecule- fermion correlation, with the associated density Pmf- 



Reference 56J outlined the limitations of the non-equilibrium theory, by claiming that to 



obtain the correct two-body physics in the low-density limit it would be necessary to include 
three point and possibly higher correlations. While this fact is indeed true, we have amended 
it in the previous section. For experimentally reasonable parameters, the mean-field theory 
can be complemented by the correct renormalized propagator to accurately describe the 



26 



equilibrium properties of the gas. Encouraged by this argument, we now apply it to the 
problem of a field ramp as well. 

In the following we wish to study molecular formation via a time dependent ramp of the 
magnetic field across the resonance. To this end we use two approaches: the first consists 
of propagating equations (l33.aH33.dl) . ramping the detuning linearly in time from a large 
positive value to a large negative one, and plotting the final molecular population as a 
function of detuning ramping rate R. The second approach consists in noticing that if u(t) 
is a linear function of time, then the mean- field Hamiltonian (eq. 1221) is ideally suited to a 
Landau-Zener treatment, whereby the final molecular population as a function of detuning 
can be readily written as 

n m /min(n b ,nf) = 1 — (34) 

Here n m /max{ni,,nf) is the fraction of possible molecules formed, and R = is the 

inverse ramp rate, and the exponential time constant is given by r = 4^- = ha ^ f b A B > wri ere 
a^g is the background scattering length, h is Plank's constant, and A B is the magnetic field 
width of the resonance. 

Remarkably, the characteristic sweep rate r does not depend on the fermionic density. 
This arises from the fact that in mean-field theory the momentum-states of the fermionic 
gas are uncoupled, except via the depletion of the condensate. Since in the Landau-Zener 
approach the depletion is assumed small, it follows that the various fermionic momentum 
states are considered independently, and thus the probability of transition of the gas is equal 
to the probability of transition of each individual momentum state. This approximation is 
only valid for narrow resonances, such as the 492. 5G resonance in table HI 

Molecular formation rate versus ramp rate is shown in Figure [121 f°r the cases rbf = 6 
(more fermions than bosons) and r&/ = 0.6 (more fermions than bosons. In both cases 
the Landau-Zener result agrees nearly perfectly with direct numerical integration. This 
agreement is surprising in the case of more fermions, since the width of the crossing is 
proportional to density of leftover bosons, and we expect that this number will change 
substantially as the bosonic population is depleted via the formation of molecules. This 
type of time dependent crossing should not be properly described by the Landau-Zener 
formula. However, monitoring the time evolution of the molecular population as a function 
of time shows that the majority of the transfer takes place quite abruptly somewhat after 
crossing the zero-detuning region, whereby the change in bosonic density does not modify 
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the energy levels substantially. 




Ramp Rate (/js/G) 




500 600 



Ramp Rate (/js/G) 



FIG. 12: Transition probability into molecular state via a magnetic field ramp across the 492. 5G 
Feshbach resonance. The solid lines are obtained via numerical solutions of equations (|33.ati33.d]) . 
while the dashed lines represent the Landau-Zener equivalent. In the top panel the gas is composed 
of more bosons than fermions rj,f = 6, while in the bottom panel the opposite is true r^f = .6 
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VI. CONCLUSION 



In this article, we developed and solved a generalized mean-field theory describing an 
ultracold atomic Bose- Fermi mixture in the presence of an interspecies Feshbach resonance. 
The theory is "generalized" in the sense that it correctly incorporates bosonic fluctuations, 
at least to the level that it reproduces the correct two-body physics in the extreme dilute 
limit. This theory, like any mean-field theory, presents undeniable limitations. Nevertheless, 
any useful many-body treatment must start from a well conceived mean-field theory. Future 
directions of this work should include the generalization to finite temperature, and the 
inclusion of a trap, initially in a local-density approximation. These advances would be 
essential to check for empirical confirmation of the theory. 
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APPENDIX A: GREEN FUNCTION METHODS FOR FERMIONS 

In this Appendix we briefly introduce some of the Green function techniques that we 
found useful in our calculations. 

1. Free Green functions 

We start from the Green function for a gas of free fermions, which is given, in the 
frequency-momentum representation by 



where £(p) = p 2 /2m — /i. The momentum distribution, at equilibrium, is given by 



Here the limit comes from the equilibrium condition; the frequency, in the Green function 
definition, is the fourier space equivalent of time, whereby the real time green function 
represents the evolution of the system from from time t to t', and the observables obtained 




(A-1) 




(A-2) 
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this way, represent expected values of the kind {ip{t)\0\ip{t')) . However, since we want 
equilibrium conditions, we need to take the limit t — > t', which is non trivial, since G° is 
defined by a green function equation of the form C G°(t — t') oc 5(t — t'), where C is some 
linear operator, and which highlights a peculiar behavior in the limit we desire. However, 
since we know on physical grounds that observables, such as the momentum distribution, 
must be defined and well behaved at equilibrium, then by first taking the expectation value 
integral, and then the limit, we can circumvent the problem. In eq. IA-21 this implies we 
cannot quite get rid of the fourier transform exponent e* 1 ^* - *') until after the to integral. 

To perform the integral in IA-21 we exploit the fourier exponent, by noting that since i] is 
positive, e luJV — > as ui — > -Moo, so that the integral is identical to a contour integral over 
the path defined by the real lu axis, closed in the upper complex uj plain by an infinite radius 
semicircle, which, as we have just seen, gives no contribution to the integral. We can now 
integrate using the residue theorem. 

We note that the integrand in IA-21 has a simple pole at u = £(p) — ir] sign(^(p)). Thus, 
if £(p) > 0, then the pole is in the lower complex plane, and the integral vanishes, and if 
£(p) > 0, then the pole is in the upper complex plane, with residue 1. Using the residue 
theorem, and summarizing these results we finally get 

n(p) = ©(-£(P)), (A-3) 
which we recognize as the zero temperature fermi distribution. 



2. Interacting Green functions 

According to Dyson's equation, the green function for an interacting system has the form 

G(w, q) = 1 -, (A-4) 

w-f(p) - S 0,P) 

where T,{u, p) is an arbitrarily complicated function summarizing all the interactions in the 
system, which is known as self energy. 

The prescription to find £ is quite straightforward, and it consists of adding all amputated 
connected feynman diagrams for the system. The fact that, in general, the number of such 
diagrams is infinite, makes this task virtually impossible. Nonetheless, eq. IA-41 is very 
powerful, since it allows one to include the effect of infinite subsets of the total number of 
diagrams in the system, by only having to explicitly calculate a few representative ones. 
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An alternative standard approach, leads to the exact result (NOTE: Abrikosov measures 
energy from /x, here we measure from 0, which is more standard.) 



G(u,p) 



duj' 



A(uj,p) B(u,p) 



(A-5) 



UJ — Uj' + IT] UJ + Uj' — it] 

Where A and B are, again, arbitrary complicated functions, though they are known to be 
finite. 

To understand A and B more closely, we need to introduce the following well known 
identity: 



lim 



1 



— = V— =F i^S(x) 
^->0 x ± iv x 



(A-6) 



where V is a Cauchy principal value, which represents the contribution due to a discontinuity 
in a Riemann sheet (branch cut), and the delta function represents the contribution due to 
the pole. 

Applying IA-61 to [A-51 we get 



POO 

Re G(u, p) = V duj' 
Jo 

Im G(u, p) = 



A(u,p) B(u,p) 



UJ — UJ + ITj UJ + UJ' — IT] 

—7tA(uj, p) if uj > 
ttB(—uj, p) if uj < 



(A-7) 
(A-8) 



Finally, eq. IA-21 represents a fundamental property of green functions, and it can be 
generalized to interacting systems simply substituting G° with G. Applying it to eq. I A- 51 
and performing the uj integral first, we get 



n(p) 



duj'B(uj, p) 



duj — Im G(uj,p). 
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(A-9) 



Introducing the function p(uj,p) = —2ImG(uj,p), generally called spectral function, the 
above equation can be written as 



n 



'JP) = J ^P(^>P)©(-^)- 



(A-10) 



An important property of the spectral function is that for all p, 

duj 



2tt 



(A-ll) 



This can be understood as a sum rule in the following sense: if we wish to calculate the 
number of holes in the system, we would take the rj —>■ 0" limit in equation IA-21 The 
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distribution would then have been n/ lo ; es (p) = 1 — n(p) = J" °° du'A(u', p), so that 1 
/ °° du/ [A(u/, P ) + B(u>, p)] = / f p(u,, p). 

Using eq JA-4l together with the definition of p, we can write 



p(u;,p) = - ( A " 12 ) 

[co> - £(p) - 7?eE(u;, p)] + [7mE(w, p)j 

Furthermore, if £ were to be real, or if, equivalently, the pole of the green function were to 
be real, for some momentum p, then taking the limit 7mS — > of IA-41 and using eq lA-6l 

we get 

p(u, p) = 2n5(co - £(p) - i?e£(cu, p)), (A-13) 
which can be simplified, using the properties of the delta function, to 

p(u, p) = 2<kZ{p)5{u - u (p)), (A-14) 

where Z, known as spectral weight is given by 

Z (P) = h 9 B M ' (A " 15) 

1 — -^--KeLiu;, p) , . 
and o; (p) is the pole of the green function, defined by 

wo(p) -e(p) - S(wb(p),p) = 0. (A-16) 

The momentum distribution in this case is thus given by 

n(p) = Z(p) j 5{uj - u {p))e{-u) = Z(p)Q(-Lo (p)). (A-17) 
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